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Neutron scattering experiments were performed to investigate magnetic order and magnetic excitations in 
ternary iron chalcogenide K2Fe4Se5. The formation of a superlattice structure below 580 K together with the 
decoupling between the Fe-vacancy order-disorder transition and the antiferromagnetic order transition appears 
to be a common feature in the A2Fe4Se5 family. The study of spin dynamics of K2Fe4Se5 reveals two distinct 
energy gaps at the magnetic Brillouin zone center, which indicates the presence of magnetic anisotropy and the 
decrease of local symmetry due to electronic and orbital anisotropy. The low-energy spin wave excitations of 
K2Fe4Se5 can be properly described by linear spin wave theory within a Heisenberg model. Compared to iron 
pnictides, K2Fe4Se5 exhibits a more two-dimensional magnetism as characterized by large differences not only 
between out-of-plane and in-plane spin wave velocities, but also between out-of-plane and in-plane exchange 
interactions. 

PACS numbers: 74.70.Xa, 75.25.-j, 75.30.Ds, 78.70.Nx 



New excitement in research on iron-based superconductors 
has arisen recently due to the discovery of the new super- 
conducting compound KxFe2-y^^2 with Tc above 30 K 1 1]. 
Superconductivity with similar critical temperature has been 
found soon after in other isostructural AxFe2-y^^2 compounds 
with A = Rb, Cs and Tl ||2|-[4|]. The extraordinary character- 
istic of electronic band structure in this system is the absence 
of hole pockets at the zone center. This poses a strong chal- 
lenge to the well accepted paradigm for Fe-based supercon- 
ductors concerning the pairing symmetry and the nature of 
magnetism based on the Fermi surface nesting scenario [5\. 
Besides, the V5 x V5 type of Fe-vacancy order was observed 
and it was accompanied by the formation of antiferromag- 
netic blocks of Fe spins mM\. The optimal composition of 
A2Fe4Se5 is suggested according to the observed ^f5x ^f5 Fe- 
vacancy-order pattern. In addition to the ^f5x ^f5 superstruc- 
ture phase, detailed experiments using transmission electron 
microscopy, x-ray, and neutron scattering revealed the exis- 
tence of a V2 X V2 superstructure phase in superconducting 
samples [8, 9]. It is suggested that superconductivity might be 
present in the ^f2x V2 phase rather than the V5x V5 phase, 
which has a Mott insulator ground state. Further experimental 
work provided evidence of the phase separation on nanoscopic 
length scales in AxFQ2-y^^2 compounds |10-12]. A thorough 
understanding of the physical properties of the superconduct- 
ing phase as well as the relation between the superconducting 
and the antferromagnetic phases is still needed. 

Although the V5 X V5 Fe-vacancy-ordered phase is a 
non- superconducting phase, it does exhibit close relationship 
with the superconducting phase. Actually, the V5 x V5 Fe- 



vacancy-ordered phase appears inevitably in all AxFQ2-y^^2 
bulk superconducting compounds. Given the fact that dy- 
namic magnetism in Fe-based superconductors may play an 
important role in mediating superconductivity [Il3l] . it is nec- 
essary to understand both static and dynamic magnetism in 
the insulating V5 x V5 Fe-vacancy-ordered phase. The spin 
wave of block antiferromagnetic A2Fe4Se5 phase has been 



examined theoretically by many groups [14l-ll7l]. Further- 
more, time-of-flight inelastic neutron scattering experiments 
have also been performed to investigate the spin wave disper- 
sion of insulating Rbo.89Fei.58Se2, (Tl,Rb)2Fe4Se5, and super- 
conducting Cso.8Fe1.9Se2 compounds [18-20]. These experi- 
ments reveal that the antiferromagnetic next-nearest-neighbor 
couplings exhibit comparable strengths in both iron pnictides 
and iron chalcogenides [21, i2] . 

In contrast to the neutron time-of-flight technique, triple- 
axis neutron spectrometry can provide more accurate dynamic 
information in a given region of energy and momentum space. 
In the present work, we studied the magnetism, especially the 
low-energy spin wave dispersion, of a block antiferromagnetic 
K2Fe4Se5 compound with both difl'raction and triple-axis in- 
elastic neutron scattering techniques. We find that the energy 
spectrum at the Brillouin zone center can be modeled with 
two diff'erent spin anisotropy gap parameters, which can be 
interpreted as indications for a breaking of local symmetry. 
A detailed analysis of the spin wave dispersion shows that 
the low-energy magnetic excitations of the antiferromagnetic 
block phase of K2Fe4Se5 can be well described by a Heisen- 
berg model with local magnetic exchange couplings extended 
to the third-nearest-neighbor. The larger energy band width 



observed in K2Fe4Se5 is related to a stronger exchange cou- 
pling strength and higher antiferromagnetic transition temper- 
ature compared to other AxFQ2-y^^2 compounds. 

The single crystals of K2Fe4Se5 were grown by the Bridg- 
man method as previously reported [23]. Single-crystal neu- 
tron diffraction and inelastic neutron scattering measurements 
were performed on the thermal-neutron two-axis diffractome- 
ter D23 and triple-axis spectrometer INS at the Institut Laue 
Langevin (Grenoble, France). The crystal we used for the neu- 
tron scattering measurements was in the shape of a cylinder 
with a total mass of 3.9 g. For the diffraction measurement 
at D23, a Cu (200) monochromator was selected to produce a 
monochromatic neutron beam with wavelength of 1 .28 A. For 
the inelastic neutron scattering measurements at INS, pyrolitic 
graphite PG(002) was chosen as analyzer, while Si(lll) or 
Cu(200) was selected as monochromator. The fixed final neu- 
tron energy of £"/ = 14.66 or 34.83 meV was used depending 
on different momentum and energy ranges. For convenience, 
in this paper we will always describe the neutron scattering 
data in the high symmetry tetragonal I4/mmm space group no- 
tation with momentum transfer wave vectors Q = (qx,qy,qz) (in 
units of A"^) at position (HKL) = (qxa/ln, qyb/ln, qzcl2n) in 
reciprocal lattice units, where a = b = 3.898 A and c = 14.121 
A at r = 2 K. The K2Fe4Se5 single crystal was aligned in the 
scattering plane defined by the orthogonal vectors (2 1 0) and 
(0 1), in which spin wave excitations along main symmetry 
directions in the magnetic Brillouin zone can be surveyed. 

The block antiferromagnetic structure [Fig. 1(a) and (b)] 
with the formation of the V5 x V5 superlattice has been sug- 
gested to describe the low-temperature structure for A2Fe4Se5 
(A = K, Rb, Cs andTl) systems L^,i7,]- Each block is composed 
of four Fe spins and these spins aligned ferromagnetically 
within the block, while the spins of neighboring blocks are 
aligned antiferromagnetically. Our neutron scattering results 
again confirmed the block antiferromagnetic structure as indi- 
cated by the contour map measured in (HKO) and (HKl) plane 
of the I4/mmm unit cell [Fig. 1(c) and (d)]. The magnetic 
structure refinement based on collected integrated intensities 
gives ordered moments of 3.2(3) jib and 2.7(3) jib at 2 and 300 
K, respectively. In Fig. 1(e), the rapid increase of the inten- 
sity of (3/5 1/5 0) structural peak below Ts = 580(3) K indi- 
cates the structural phase transition from the high- temperature 
Fe- vacancy-disordered phase with I4/mmm symmetry into the 
low-temperature Fe- vacancy-ordered phase with I4/m symme- 
try. The temperature variation of the integrated intensity of 
(2/5 1/5 1) magnetic peak can be fit with an empirical power 
law / oc (1-T/TNf^, which yields Tn = 553(3) K and a critical 
exponent yS of 0.27(1), as represented by the solid line in Fig.l 
(e). 

In Fig. 2(a), the energy scan at the Brillouin zone center Q 
= (2/5 1/5 -3) demonstrates a clear energy gap of magnetic ex- 
citation due to single ion anisotropy. The significant spectral 
signal is observed to extend up to 20 meV after substraction 
of background that was obtained by a comparable scan at Q = 
(0.7 0.35 -3). The existence of an energy gap is further demon - 
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FIG. 1: (Color online) (a) Crystal and magnetic structure of 
K2Fe4Se5 insulating phase. Gray lines highlight the low symmetry 
tetragonal unit cell, (b) Vertical view of antiferromagnetic spin block 
in Fe layer. The Vs x Vs superlattice unit cell is marked by gray 
bonds, while the dotted lines mark the high- temperature I4/mmm unit 
cell. The blue and yellow spheres denote Fe atoms with spin-up and 
spin-down moments. The exchange couplings (/i, J2, J3, /J, J2 and 
73) indicate intra- and inter-block exchange interactions, (c) and (d) 
Experimental contour map in the first quadrant of (HKO) and (HKl) 
plane in reciprocal space, where nuclear reflections and magnetic re- 
flections from both chiral domains are observed, (e) Temperature 
dependence of integrated intensities of (3/5 1/5 0) nuclear and (2/5 
1/5 1) magnetic superstructure reflections. 



strated by Q- scans along two principal momentum directions 
below and above the gap energy as shown in Fig. SI (a) and 
Fig. S2(a). To determine the accurate gap energy, the ob- 
served energy scan data was used to fit within the framework 
of empirical spin wave relation convoluted with instrumen- 
tal resolution function |l24 U5] . Interestingly, the data in the 
whole investigated energy range cannot be properly fitted with 
only one energy gap parameter. A second energy gap param- 
eter has to be introduced. As a result, the analysis leads to 
two gap energies of 8.7(3) and 16.5(3) meV. A similar feature 
was recently reported for BaFe2As2 where inelastic polarized 
neutron scattering studies lead to the observation of a strongly 




Energy (meV) 




17- 


r\ 4 


L [ ' 


^ 16- 
> 

CD 

E 15- 


t t i^" 


n 






■^1 


> J 

CD 9- 
LU 


-»- Energy Gap 2 
-•-Energy Gap 1 


' 


H^i 


I 




8- 


^r^ 


H 




^ 1 


7- 




' 1 



4 8 12 

Energy (meV) 



100 200 

Temperature (K) 



FIG. 2: (Color online) (a) Energy scan at the antiferromagnetic wave 
vector Q = (2/5 1/5 -3) at 2 K. The solid line is the fitting result 
with the dispersion model convoluted with the instrument resolution 
function. The spectral weights arising from two energy gap terms 
are highlighted as the shaded areas under the curve, (b) Energy scan 
at the zone boundary with wave vector (2/5 1/5 -4) at 2 K. (c) Tem- 
perature dependence of energy scan at wave vector Q = (2/5 1/5 -3). 
Solid lines indicate fitting results where corrected Bose population 
factor is also taken into account. Background scattering data (BKG) 
at an arbitrary wave vector is also shown for comparison. The curves 
are shifted upward to clearly show the change in energy gap features, 
(d) Variation of the two gap energies in the investigated temperature 
range. 
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FIG. 3: (Color online) (a) Spin wave dispersion relation along L di- 
rection as deduced from Q- scans at constant energy. The solid line 
is the fitted dispersion with a Heisenberg model as described in the 
text, (b) Variation of dynamic spin structure factors, (c) Spin wave 
dispersion along (2// H 0) direction. Solid lines are the fitted spin 
wave dispersions originating from two diff'erent structural domains 
(Dl and D2) with equal populations, (d) Variation of dynamic spin 
structure factors along (2// H 0) momentum space direction around 
Q = (2/5 1/5 1). 



anisotropy spin excitation and the occurrence of two gap ener- 
gies 1 26] . The observed strong in-plane single ion anisotropy 
in BaFe2As2 demonstrated the important role of orbital de- 
grees of freedom in the iron pnictides. Regarding A2Fe4Se5, 
it was found that the magnetic exchange energy will be min- 
imized and the electron correlation will be enhanced by the 
presence of vacancy order in the V5 x V5 superlattice phase. 
Furthermore, a particular orbital order pattern will be pre sent 
and it will result in the breaking of local fourfold symmetry 
on the Fe site ll27n . The observed two energy gap features 
in K2Fe4Se5 might be an indication of the reduction of local 
symmetry due to electronic and orbital anisotropy. Despite the 
two energy gaps, only the one dominating the spectrum was 
used for the subsequent spin wave analysis. It is also notewor- 
thy that the gap energy values obtained in K2Fe4Se5 are of the 
same order of magnitude as the ones of various iron arsenides, 
in spite of quite large differences in ordering temperature and 
magnitude of Fe moments [28]. 

Similar fitting procedures are performed to fit energy scan 
data observed at different temperatures [Fig. 2(c)]. The ob- 
tained two spin gap values as a function of temperature are 
shown in Fig. 2(d). Both energy gaps exhibit similar temper- 
ature dependence in spite of their diflTerent spectral weights. 
Narrowing of the energy gaps, which demonstrates the reduc- 



tion of spin anisotropy, is clearly observed upon the increase 
of temperature from 2 to 300 K. Apart from the spin gap at 
the zone center, the Brillouin zone boundary is also examined. 
Fig. 2 (b) shows an energy scan at Q = (2/5 1/5 -4) located 
at the magnetic Brillouin zone boundary. The zone boundary 
energy is deduced to be 25.3(3) meV and is almost three times 
larger than the gap value at zone center. Spin wave energies at 
both zone center and zone boundary will help us to model the 
spin wave dispersion relationship in K2Fe4Se5 as discussed in 
the following text. 

To determine the spin wave dispersion in the block antifer- 
romagnetic K2Fe4Se5 phase, constant energy measurements 
around the zone center Q = (2/5 1/5 1) were performed along 
two high symmetry directions in the Brillouin zone, i.e. (0 
L) and (2// H 0) directions (Fig. SI and Fig. S2). The disper- 
sion relations extracted from the measured scans are plotted 
as spheres in Fig. 3 (a) and Fig. 3 (c). Given the fact that 
K2Fe4Se5 is an antiferromagnetic insulator with a rather large 
ordered Fe moment, we analyze spin wave dispersion data in 
linear spin wave approximation with Heisenberg Hamiltonian 
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FIG. 4: (Color online) (a) and (b) Low energy spin wave branches 
along the selected high symmetry directions in the magnetic Bril- 
louin zone. Dl and D2 represent spin waves originated from two dif- 
ferent chiral domains, (c) Illustration of the defined high- symmetry 
points in the Brillouin zone. P = (0, 0, In/c); M' = (4;7r/5a, In/Sa, 
In/c); X' = (3;r/5a, -n/Sa, In/c); Y' = (n/Sa, 3;r/5a, In/c); T! = 
(0, 0, 4;r/c); A' = {Anj'Sa, In 1 5a, Anic). (d) Calculated low energy 
spin wave density-of- states for isostructural compounds K2Fe4Se5 
andRbo.89Fei.58Se2. 



coupling constants (/c), while Ds is the uniaxial single-ion 
anisotropy constant. 

As illustrated in Fig. 1(b), eight Fe spins are distributed 
in the magnetic unit cell of I4/m symmetry, thus one two- 
fold degenerate acoustic mode at lower energy together with 
three two-fold degenerate optical modes at higher energy are 
expected. As observed in isostructural Rbo.89Fei.58Se2 and 
(Tl,Rb)2Fe4Se5 by time-of-flight neutron scattering measure- 
ments, all three optical modes are located at energies higher 
than 100 meV fis', 'l9']. In our work, a gapped acoustic 
mode is measured accurately, but unfortunately the optical 
modes could not be reached because of the kinematic limit 
and low neutron flux at high transfer energy. The acoustic 
mode is found mainly due to the antiferromagnetic interaction 
between ferromagnetic blocked spins and it is mainly deter- 
mined by J[, J^, Ji,, Jc and Dg within the J1-J2-J3 model. For 
instance, the energy gap at the zone center can be expressed as 
As = S ^Ds(2J[ -\- 4/2 + 4/3 + Ds), while the spin wave en- 
ergies at the zone boundary Q = (2/5 1/5 2) and at Q = (0 1) 
are found to be Ezb = S Y[2(7[T27^T273yTDn(477TD;) . 
Therefore, we adopted the exchange parameters S Ji = -36(2) 
and SJ2 = 12(2) meV from isostructural Rbo.89Fei.58Se2 lilSll 
and fit our observed spin wave dispersion data with /J, J2, 
/s, Jc and Ds as variables. The spectra originated from two 
diflferent chiral domains were also taken into account. The 
fitting results exhibit reasonable agreement with experimental 
data and yield exchange parameters as S J[ = 11(2) , S J!^ = 
19(2), SJ3 = 12(2) ,SJc = 0.88(8) and ^Z), = 0.46(6) meV. 



In the V5 X V5 vacancy ordered A2Fe4Se5 phase, the near- 
est and the next-nearest neighbor interactions (/i, /J, J2, J 2) 
are undoubtedly significant to describe its spin wave behav- 
iors. However, the importance of the third-nearest-neighbor 
interaction (J3) is still a matter of controversy U^ 19]. In 
the present work, we found that the influence of J3 reflected 
mainly in the change of the zone boundary energy of the 
acoustic branch. Based on the analysis on our own data 
as well as the results presented in Ref. |[18l ll9|], we argue 
that the third-nearest-neighbor exchange interaction J3 also 
plays a role in determination of spin dynamics in A2Fe4Se5. 
By considering the interactions extended to the third-nearest- 
neighbor, variations of the dynamic spin structure factor 
S(Q,<jj) are obtained and presented in Fig. 3(b) and (d) to 
illustrate the intensity distribution and the individual contri- 
bution from two diflTerent domains. It is clearly shown that the 
dynamic structure factors from two domains exhibit closeness 
and even overlap in certain Q range in the (2H H 0) direction. 

The spin wave dispersion of the acoustic modes along all 
high symmetry directions is calculated by using the obtained 
exchange parameters, as shown in Fig. 4(a) and (b). The spin 
wave density-of- states (SWDOS) of K2Fe4Se5 was also ob- 
tained by the summation over all wave vectors in the Brillouin 
zone, as plotted in Fig. 4(d). The density-of- states directly 
reflect the distribution of spin wave energies. Similar to pnic- 
tides, SWDOS of K2Fe4Se5 exhibits sharp peaks or anomalies 
at energies of van Hove singularities [29, 30]. The SWDOS 
of Rbo.89Fei.58Se2, calculated using the exchange parameters 
provided in Ref. [18j, is also given in Fig. 4(d) for compari- 
son. It can be seen that the energy bandwidth of K2Fe4Se5 is 
larger than that of Rbo.89Fei.58Se2, which reflects stronger ex- 
change interactions and a higher antiferromagnetic transition 
temperature. 

The investigations of the low energy spin wave disper- 
sions of the parent phase of iron pnictides, e.g. in BaFe2 As2, 
CaFe2As2 and SrFe2As2, have revealed values of vjva = 0.2- 
0.5 for the ratio between out-of-plane and in-plane spin wave 
velocities (vJVa) [131 l25l |3l|]. Therefore, it was suggested 
that the parent phases of iron pnictides are three dimensional 
antiferromagnets rather than quasi two-dimensional antifer- 
romagnets as layered cuprates, in which finite but small in- 
terlay er coupling exists [32]. The spin wave dispersion of 
K2Fe4Se5 at lower energy yield an in-plane velocity of Va = 
380(20) meV-A and an out-of-plane velocity of v^ = 50(10) 
meV-A. The ratio of vJVa = 0.13(2) observed in K2Fe4Se5 
indicates that it possesses more two-dimensional character in 
magnetism than iron pnictides. The more two-dimensional 
magnetism in K2Fe4Se5 is also reflected by the ratio between 
out-of-plane and in-plane exchange interactions. For instance 
a ^"^-"^-^^^"^ ratio of less than 3% is obtained if we compare 

J in -plane 

Jc with nearest neighbor exchange interaction /i or eflTective 
block in-plane exchange parameter Jeff = \(J[-^ 2/2 + 2/3) as 
suggested in Ref. [18]. It is believed that superconductivity is 
favored in more two-dimensional magnetic system since large 
spin fluctuations might suppress the long range order and me- 



diate superconductivity. 

In summary, we have investigated both static and low 
energy dynamic magnetism of a block antiferromagnetic 
K2Fe4Se5 single crystal by using elastic and inelastic neu- 
tron scattering techniques. K2Fe4Se5 is found to undergo 
a Fe- vacancy order-disorder transition and an antiferromag- 
netic transition at Ts = 580(3) and Tn = 553(3) K, respec- 
tively. Additionally, we obtained the acoustic mode of spin 
wave dispersions and it was properly fitted by using a linear 
spin wave model based on J1-J2-J3 Heisenberg exchange cou- 
plings. Moreover, we found two well separated energy gaps 
at the magnetic Brillouin zone center. The breakdown of lo- 
cal tetragonal symmetry due to the emergence of orbital order 
might be responsible for the appearance of two energy gaps. 

T.W. thanks the Deutsche Forschungsgemeinschaft for fi- 
nancial support under the DFG Priority Program 1458. 
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